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> ABSTRACT 

o 

Context. Magnetic clouds (MCs) are a subset of interplanetary coronal mass ejections (ICMEs) where a magnetic flux rope is detected. 
Is the difference between MCs and ICMEs without detected flux rope intrinsic or rather due to an observational bias? 
Aims. As the spacecraft has no relationship with the MC trajectory, the frequency distribution of MCs versus the spacecraft distance 
to the MCs axis is expected to be approximately flat. However, Lepping and Wu (2010) confirmed that it is a strongly decreasing 
function of the estimated impact parameter. Is a flux rope more frequently undetected for larger impact parameter? 
Methods. In order to answer the questions above, we explore the parameter space of flux rope models, especially the aspect ratio, 
boundary shape, and current distribution. The proposed models are analyzed as MCs by fitting a circular linear force-free field to the 
magnetic field computed along simulated crossings. 
Results. We find that the distribution of the twist within the flux rope, the non-detection due to too low field rotation angle or 
magnitude are only weakly affecting the expected frequency distribution of MCs versus impact parameter. However, the estimated 
Oh impact parameter is increasingly biased to lower values as the flux-rope cross section is more elongated orthogonally to the crossing 
JL trajectory. The observed distribution of MCs is a natural consequence of a flux-rope cross section flattened in average by a factor 

2 to 3 depending on the magnetic twist profile. However, the faster MCs at 1 AU, with V > 550 km/s, present an almost uniform 
distribution of MCs vs. impact parameter, which is consistent with round shaped flux ropes, in contrast with the slower ones. 
Conclusions. We conclude that either most of the non-MC ICMEs are encountered outside their flux rope or near the leg region, or 
they do not contain any. 
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> 
CO 

^j" 1. Introduction ping et al. 1990; Lynch et al. 2003; Dasso et al. 2006; Leitner 

CO et al. 2007). 

Interplanetary Coronal Mass Ejections (ICMEs) are detected in 
^— H the solar wind (SW) by in situ plasma and magnetic field mea- 

t-H surements onboard spacecraft. They are the counterpart of Coro- Gosling (1990) found that a MC is present inside ICMEs 

CS| nal Mass Ejections (CMEs) observed with coronagraphs (e.g. only for 10-30% of the cases. Presently, in average, a MC is de- 

Howard 201 1 ; Lugaz & Roussev 201 1). With STEREO twin tected in about 30% of ICMEs (Richardson & Cane 2010; Wu & 

spacecraft having both in situ and imager instruments, this link is Lepping 201 1). Cane & Richardson (2003) found that this ratio 

• i-H presently well etablished (e.g. Harrison et al. 2009; Kilpua et al. is evolving with the solar cycle: the MC/ICME ratio increases 

^2011; Rouillard 2011; Lugaz et al. 2012; Wood et al. 2012, and from « 15% at solar maximum to « 100% at solar minimum. 

H references therein). ICMEs are defined by one or several criteria They interpreted this evolution as due to an observational selec- 

^ (see Wimmer-Schweingruber et al. 2006; Zurbuchen & Richard- tion effect since CMEs are launched from higher solar latitudes 

son 2006, for reviews). Typical criteria are: (a) a stronger mag- at solar maximum than at minimum, so a spacecraft located in 

netic field with lower variance than in the surrounding SW; (b) the ecliptic more frequently crosses the flux rope away from its 

a low proton plasma p p (< 0.4 typically); (c) a smooth and large axis at solar maximum than at solar minimum (Richardson & 

rotation of the magnetic field; (d) a proton temperature at least Cane 2004). This evolution with solar cycle is confirmed by 

lower by a factor 2 than in ambient SW with the same velocity as their newer results (Richardson & Cane 2010) and by the results 

in the MC; (e) an enhanced helium abundance (He/H > 6%); (f) of Kilpua et al. (2012) around the minimum between solar cycle 

the presence of counter- streaming suprathermal (> 80 eV) elec- 23/24 (they found that « 76% of ICMEs have flux rope charac- 

tron beams; (g) enhanced ion charge states. Magnetic clouds teristics in the time period 2007-2010). Since nearly all ICMEs 

(MCs) are defined with criteria (a-d) all satisfied (Burlaga et al. are MCs at solar minimum, it has been suggested that MCs are 

1981), then MCs are a sub-class of ICMEs. This signature has only observed when the spacecraft crosses the magnetic struc- 

been modeled as a magnetic flux rope (e.g. Burlaga 1988; Lep- ture near the flux rope center (e.g. Jian et al. 2006). 
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Fig. 1. Probability distribution, P bs(p), of the impact parameter 
(p). The results of MCs observed by WIND at 1 AU and fitted by the 
Lundquist model (Lepping et al. 1990; Lepping & Wu 2010), are shown 
with a histogram having 10 bins of p. Black curve: A linear fit to the 
histogram. Red curve: The gaussian function derived by Lepping & Wu 
(2010). The histogram in the top panel has 100 MCs, while the one in 
the bottom panel is restricted to the 74 best observed MCs (quality 1,2). 



Still, the in situ observations provide only a ID cut through 
a 3D structure so there is a lack of information. For MCs this is 
typically complemented by a fit of a magnetic model to the data, 
providing both the local orientation of the flux rope and its field 
distribution within the cross-section. The most often used model, 
so far, is the so-called Lundquist model (Lundquist 1950), which 
considers a static and axi-symmetric linear force-free magnetic 
equilibrium configuration (e.g. Goldstein 1983; Burlaga 1988). 
Its main advantage is its simplicity (low number of free parame- 
ters) while it satisfies the low plasma f3 condition typically found 
in MCs (typically j3 p < 0.1) and fits relatively well observations 
(e.g. Lepping et al. 1990; Burlaga 1995; Lepping et al. 2003; 
Das so et al. 2005b). The fit of the Lundquist model provides an 
estimation of the closest approach position (CA) of the space- 
craft trajectory to the flux rope axis. It is generally expressed in 
% of the flux rope radius R (e.g., Lepping & Wu 2010). CA/R is 
also called the impact parameter, and noted p (e.g. Lynch et al. 
2003; Jian et al. 2006). The sign of p indicates which side of the 
MC is crossed by the spacecraft. Below, we consider only the 
distance to axis, so \p\, and we simplify the notation to p. 

With a set of 98 MCs observed at 1 AU, Lepping & Wu 
(2010) found that the number of MCs, detected at 1 AU nearby 
Earth, decreases with |CA|, or p (Fig. 1, Sect. 2.1), confirming 
previous results (Lepping et al. 2006). They checked that the de- 
pendence with p of the rotation angle and of the mean magnetic 
field along the spacecraft trajectory was behaving as expected for 
the Lundquist field model. While p is the most uncertain param- 
eter of the fit result for an individual MC (Lepping et al. 1990), 
this self-consistency tests reinforce that, in average, p was esti- 



mated correctly enough by the fit of the Lundquist model to the 
observations. 

Every MC with a flux-rope axis inclined on the ecliptic plane 
crosses it. Then, because CMEs depart from the Sun at any lon- 
gitude relative to the center disk, the related MC in situ observa- 
tions are expected to cross the flux rope at random distance from 
its axis. Next, we consider the minority of MC cases where the 
flux-rope axis is nearly parallel to the ecliptic plane. Because 
CMEs depart typically away from the solar equator, the space- 
craft is expected to cross the flux rope at a distance to its axis 
which is correlated to its launched latitude. Such cases imply a 
bias towards a larger impact parameter (even if the deflection 
of CMEs toward the heliospheric current sheet decreases this 
effect). From these considerations, one expects a flat, or even 
slightly increasing, distribution of MCs versus p which is not 
observed (Fig. 1). 

A first interpretation of the observed decrease (Fig. 1) is a 
strong selection effect due to a greater difficulty to detect a flux 
rope when p is larger. In this case, correcting this selection ef- 
fect by supposing a flat distribution, with the frequency detected 
for low p value, would typically double the number of detected 
MCs. Then, does a large fraction of the non-MC ICMEs cor- 
respond to undetected flux rope with the spacecraft trajectory 
being too far from the flux rope axis? 

A weakness of the above analysis is that the deduced im- 
pact parameter can still be biased by the selection of a particular 
model. Indeed, the self-consistency tests of Lepping et al. (1990) 
only check that the fit to the data provides coherent results with 
the hypothesis of the model. Evidences of compression in the 
direction of propagation are present in CME observations (e.g. 
Savani et al. 2009, 2010) and in MHD simulations (e.g. Cargill 
& Schmidt 2002; Odstrcil et al. 2004; Manchester et al. 2004; 
Lugaz et al. 2005b; Xiong et al. 2006). Such a compression 
flattens the cross section, and such geometrical feature has been 
partly taken into account by Vandas & Romashets (2003). They 
developed an extension of the Lundquist model from a circular 
to an elliptical boundary. This model is still analytical (but rel- 
atively complex) and it introduces only one more parameter, the 
aspect ratio of the ellipse, if one supposes that the major axis 
of the elliptical cross section of the flux rope is perpendicular to 
the direction of its motion. Moreover, it provides a better fit to 
observed MCs having a relatively uniform field strength. This 
flatness of the magnetic profile increases with the aspect ratio, 
pointing that some MCs have a relatively flat cross section (Van- 
das et al. 2005; Antoniadou et al. 2008). Finally, Vandas et al. 
(2010) tested this model with the results of an MHD simulation 
by exploring several spacecraft crossings of the simulated flux 
rope. They concluded that both the aspect ratio and the impact 
parameter p are fully reliable only for low p values, confirming 
and extending Lepping et al. (1990) results. 

A variety of alternative models have been proposed for MCs. 
Keeping the cylindrical symmetry, a variety of non-linear force- 
free field models are also possible (see Sect. 3.1). One possibility 
has a uniform twist within the cross-section (e.g. Farrugia et al. 
1999; Dasso et al. 2003, 2005a). Also, several non force-free 
models have been applied, using different shapes for their cross 
sections (e.g. Mulligan et al. 1999; Cid et al. 2002; Hidalgo et al. 
2002; Hidalgo 2011). So far, even if a given model has been 
shown to better fit the data of a few MCs than other models, this 
conclusion has not been extended to a large set of MCs. Indeed, 
the typical internal structure (e.g. the twist profile) of MCs is 
still not precisely known. 

Another approach is to include the curvature of the flux 
rope axis by developing toroidal models (Marubashi 1997; Ro- 
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mashets & Vandas 2003; Marubashi & Lepping 2007; Ro- 
mashets & Vandas 2009; Owens et al. 2012). This approach 
is especially needed when a leg of the flux rope is crossed (i.e. 
when the spacecraft trajectory is close to the local flux-rope axis 
direction). In this case, the inclusion of the axis curvature can 
strongly change both the deduced axis orientation and impact 
parameter (Marubashi et al. 2012). Such leg crossings are typi- 
cally characterized by a long duration MC with a complex rota- 
tion profile of the magnetic field as well as a low angle between 
the solar radial direction and the flux rope axis (called cone an- 
gle). The frequency of such cases is small in the data set of Lep- 
ping & Wu (2010) with, for example only 6 over 98 MCs with 
a cone angle below 30°, and none for the 67 MCs of quality 1,2 
(as defined in their paper). Then, we consider only local models 
of MCs with straight axis as they have less free parameters than 
toroidal models. 

One major unknown is the extension of the MC cross section. 
A way to deduce it is to solve the non-linear force free equations 
by a direct numerical integration with the measured vector mag- 
netic field as boundary conditions. This approach only supposes 
a magnetostatic field invariant by translation along the straight 
axis (Hu & Sonnerup 2002; Sonnerup et al. 2006). The method 
was tested successfully with MCs crossed by two spacecraft (Liu 
et al. 2008; Kilpua et al. 2009; Mostl et al. 2009a). The results 
depend on the MC studied, ranging from nearly round to elon- 
gated cross sections (Hu et al. 2005; Liu et al. 2008; Mostl et al. 
2009a,b; Isavnin et al. 2011; Farrugia et al. 2011). The main 
limitation of this kind of approach is that it solves an ill-posed 
problem: the integration of an elliptic partial differential equa- 
tion from a part of the boundary of the domain. The results can 
indeed be strongly affected by the time resolution and range of 
the data used, as well as by the method used to stabilize the inte- 
gration. 

The short review above shows a large variety of flux-rope 
models. If the MC data of Lepping & Wu (2010) would have 
been fitted by one of the above flux rope models, how would the 
estimated impact parameter, p, been affected? Said differently, 
how model dependent is the MC distribution shown in Fig. 1? 
Moreover, how strongly does selection effects, e.g. on the 
amount of magnetic field rotation or field strength, affect such 
a distribution? We analyze these issues for a large set of MCs by 
studying a variety of force-free field models (as the plasma-/? is 
typically around 0.1 in MCs). The meaning of the main parame- 
ters used throughout the paper are summarized in Table 1 . 

The observation results and the fitting method are summa- 
rized in Sect. 2. In Sect. 3, we investigate the effect of a broad 
range of magnetic field profiles ranging from flat to peaked 
around the axis, keeping a circular cross section. In Sect. 4 
we mostly investigate the effect of the cross section elongation 
with models having elliptical cross sections. We also analyze 
the effect of bending of the cross section to a "bean shape". We 
conclude that p is most affected by the aspect ratio of the cross 
section. Then, in Sect. 5, we deduce a distribution of the aspect 
ratio compatible with the results of Lepping & Wu (2010). Next, 
in Sect. 6, we further analyze the MCs according to their global 
properties and show that some set of MCs have a relatively round 
cross-section. Finally, we conclude by summarizing our results 
and, in particular, answering the question set in the title of this 
paper (Sect. 7). 




Fig. 2. Drawing defining the geometry parameters for a spacecraft 
crossing a MC. The fit of Lundquist field is schematized by the blue 
circle while the black ellipse delimitates the half extension of the MC 
boundary. In this figure, we scale the drawing with the semi-minor axis 
of the ellipse set to unity. The true impact parameter, y/b, is larger than 
P- 

2. Observations and fitting method 

2.1. Observed probability distribution of the impact 
parameter 

We use the results of the Lundquist model fitted to MCs ob- 
served at 1 AU by WIND spacecraft from February 1995 
to November 2007. They are available in Table 2 at 
http://wind.nasa.gov/mfi/mag_cloud_Sl.html. The list contains 
the results of 120 MCs at the date of 13 Dec, 2011. However 
removing the cases where the handedness could not be deter- 
mined (flag f in the list) or the fitting convergence could not be 
achieved (flag F), this list restricts to 110 MCs. Next, we ex- 
amine the cone angle J3 which is the angle between the MC axis 
(found by the Lundquist fit) to the solar radial direction (-X axis 
in GSE coordinates). We consider a folded angle, so ft is in the 
range [0°, 90°]. Since the data obtained in the cases of MC leg 
crossing are the most difficult to analyze, making the fit results 
from cases with small ft angle the most uncertain (Sect. 1), we 
limit the study to p > 30°. This restricts the MC sample to 103 
MCs. Finally, there are 3 MCs with an impact parameter p > 1 
(so a fitted flux rope extending beyond the first zero of the axial 
field in the Lundquist model). Removing these suspicious cases, 
all of the worse class (quality 3, as defined in Lepping & Wu 
2010), it remains 100 MCs. One can even be more strict on the 
selection criteria. An extreme case is to select only the best cases 
(quality 1), with the limitation that the statistics is restricted to 
19 MCs. A less extreme case is to select the best and good cases 
(quality 1 and 2) so 67 MCs. We verify that our results are not 
significantly affected by the group of MC selected. 

Lepping & Wu (2010) found that the number of detected 
MCs decreases rapidly with p. The same result is shown in 
Fig. 1 for the 100 MCs selected and a bin size of Ap - 0.1. 
Very close results are obtained if we restrict the analysis to the 
best observed MCs (quality 1 and 2). We define a probability of 
detection by normalizing the sum of the bin counts to 1 . This al- 
lows the comparison between the results obtained with different 
sets of MCs and the model predictions. The Gaussian function 
shown in Fig. 1 (red curve) is a fit of the distribution as given 



Article number, page 3 of 15 



by Eq. (Al) of Lepping & Wu with b = and cr = 0.407. It 
also fits well both sets of MCs shown in Fig. 1. The observed 
distribution can also be fitted with a linear function without sig- 
nificant difference with a Gaussian function (taking into account 
the statistical fluctuations). 



3. Detecting circular flux ropes 

In this section, we analyze a series of circular force-free fields 
in order to test if the impact parameter could be biased by the 
choice of the model with the classical analysis of Lepping et al. 
(1990). 



2.2. Flux rope fit with the Lunquist field 

Models are used to simulate flux rope crossings, providing syn- 
thetic observations which are analyzed as MCs, so following the 
classical procedure of Lepping et al. (1990). Those synthetic 
models can be chosen as circular or elliptic, and the half exten- 
sion of such a structure is given in black in Fig. 2. Then, the 
bias in the fits are analyzed. The simulated trajectory is set at a 
distance y parallel to the x-axis (because of the invariance in the 
z direction of the models, the same B would be obtained along a 
trajectory inclined on the flux rope axis). The true impact param- 
eter is y/b where b is the size of the structure in the y-direction 
(Fig. 2). The synthetic observations are fitted with a linear force- 
free model, the classical Lundquist solution (Lundquist 1950), 
which is in cylindrical coordinates: 



B L = £ (0, Ji(ar),J (ar)), 



(1) 



where a is associated with the first zero of Bz, Bq is the axial 
field strength and J m is the ordinary Bessel function of order m. 
The fit of Eq. (1) to the synthetic observations provides an es- 
timation of y, called vl, with an origin not necessarily located 
on the true flux rope axis (Fig. 2). Then, the Lundquist fit pro- 
vides the estimated impact parameter p = y^/R where R is the 
estimated flux rope radius (defined for B z = 0). 

The B L field is fitted to the synthetic observations B b s by 
minimizing the function dev defined by: 



dev 



1 7Vp 



• Bobs,/) 2 



(2) 



where N p is the number of points in the synthetic observations, 
and B b s is related with y/b, while Bl is related to p. Providing 
that N p is large enough (i.e. N p > 20), the results of the fits 
are insensitive to the value of N p , as expected since the synthetic 
observations are well resolved with such N p values. Since the 
orientation of B b s in MCs is following that of Bl better than 
the magnetic field magnitude, Lepping et al. (1990) fitted B b s 
with B L with a two steps procedure. In a first step, both B b s and 
Bl norms are normalized to unit at each point before minimizing 
dev. Then, in a second step, the full fields are considered and dev 
is minimized by only changing the axial field strength Bq. From 
synthetic Lundquist fields, Gulisano et al. (2007) also concluded 
that fitting to normalized B b s gives better estimation of the real 
orientation of the MC axis. 

Compared to real MC observations, the models provide syn- 
thetic observations with no internal structures, and with known 
axis orientation and boundaries. The exploration of the effect 
of perturbations to B, various axis orientation and boundaries 
could be realized in the line of the following exploration of the 
parameter space (e.g. , the parameters defining the shape of the 
cross section). However, we choose to rather limit the explo- 
ration to the global structure of the flux ropes (i.e. the magnetic 
field repartition and the cross section shape) as such structure is 
expected to have a main effect on the estimated impact parame- 
ter. Then, in the above first step, the minimization is realized by 
changing p and a (as both the axis orientation and the flux rope 
boundaries, set at B z = 0, are known and fixed). 



3.1. Force-free models 

Frequently, the magnetic structure of a MC is locally approxi- 
mated by a straight flux rope invariant along its axis (Sect. 1). 
We use below an orthogonal frame, called the MC frame, with 
coordinates (x,y, z). The z direction is along the MC axis. B be- 
ing independent of z and V • B = imply together that one can 
write the magnetic field components orthogonal to the symme- 
try axis as: B x = dA/dy and B y = -dA/dx, where A(x,y) is the 
magnetic flux function. The force-free field condition implies 

aA + 



dA 



, with B Z (A) . 



(3) 



A series of non-linear force-free fields are generated by 

B z (A) = cA\ (4) 

where c and n > are independent of x, y, z. Typically, the flux 
rope boundary is set at a location where B z = which can be 
set for A = without loosing generality. We normalize the 
cross section extension to the half of its maximal value in the x- 
direction (Fig. 2). Its half maximal extension in the y-direction is 
then the aspect ratio, b, set to b = 1 in this section (i.e., circular 
shape). The flux rope axial field is called # a xis- 

Since we consider circular flux ropes in this section, Eq. (3) 
reduces to a differential equation of second order with the ra- 
dius ( ^x 2 + y 2 ). It is solved by a numerical integration using a 
shooting method (e.g. Press et al. 1992, p.746) applied to the 
resonance problem set by Eqs. (3, 4) and the three boundary con- 
ditions: A(0) = 1, [dA/dr](0) = and A(l) = (corresponding 
respectively to an azimuthal flux normalized to 1 , a regular field 
on the axis, and to B z = at the boundary). We select the lowest 
eigenvalue c to have models with an axial field vanishing only at 
the boundary as present in most MCs. For n = 1, the field is a 
linear force-free field and c - a (as defined by Eq. (1)). Finally 
the field strength on the axis can be scaled to any desired # axis 
value. 

The axial electric current density is: 

k = -A^M) = n c 2 (A(x,y)f n - 1 /^ , (5) 

where jio is the permeability of the free space. For n = 0.5, j z 
is uniform, while for n > 0.5, j z decreases from the axis to the 
boundary (where A = 0, so j z = 0). Finally, for n < 0.5, j z is 
singular at the boundary (presence of a current sheet). 

The value of n also determines the spatial variation of B. 
This is illustrated in the left panels of Fig. 3. As n is increased, 
the magnetic field strength becomes more concentrated around 
the axis and, near the boundary for n > 0.5, the azimuthal 
field (= \B y \ at y = 0) is a decreasing function of the radius 

(= ^jx 2 + y 2 ) over a larger radius range. The case n = 2, Fig. 3e, 
is an extreme case for a MC. At the opposite, as n decreases the 
profile of the magnetic field strength flattens. For n = 0.5, the 
azimuthal field is linear with radius, while for n < 0.5 it is in- 
creasing more sharply near the boundary as n is decreased. The 
case n = 0.1 is another extreme case for a MC. The Lundquist 
model fits relatively well the different models except in extreme 
cases (e.g. n = 2 case, Fig. 3e,f), even for large impact parame- 
ters (e.g. see the case y/b = 0.9 in Fig. 3b,d). 
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(a) n=0.1 , y/b=0 (b) n=0.1 , y/b=0.9 




Fig. 3. Examples of circular models (black dots) least square fitted with the Lundquist field (red curves). Three non-linear force-free models 
(n = 0.1,0.5,2) are selected to represent strong departure to the Lundquist field (n = 1). The true impact parameter, y/b, is either null (left) or 
large (right). B z is the axial field component, and B y is the azimuthal field component for y = 0. The field strength of the model on the axis, B axis , 
is normalized to 1 . 




0.2 0.4 0.6 0.8 _ 1.0 0.2 0.4 0.6 0.8 „ 1.0 




0.2 0.4 0.6 0.8 n 1.0 0.2 0.4 0.6 0.8 n 1 .0 



Fig. 4. (a) Dependence of the true impact parameter, y/b, (b) the fitted 
Lundquist field strength on the axis, B , (c) the normalized deviation, 
dev/<B>, and (d) the mean field magnitude <B> normalized to B 
in function of the impact parameter, p, found by fitting the Lundquist 
field to the models. The parameter n describes the profile of the axial 
component of the magnetic field and electric current, see Eqs. (4, 5). 

3.2. Informations provided by the Lundquist fit 

The results of Fig. 4a show that y/b is relatively well estimated 
by p for all the range of n values relevant to MCs. The extreme 
case n = 0.1 has very similar results to the case n = 0.25, so it 
is not shown. All cases with n < 1 have p only slightly lower 
than y/b so that the observed distribution probability (in function 
of p) would only be slightly compressed towards lower p values 
compared to the original distribution probability (in function of 



y/b). On the contrary, all cases with n > 1 would introduce a bias 
opposite to those observed since p > y/b (Fig. 1). We conclude 
that the deviation around y/b - p (green curve, n = 1 in Fig. 4a) 
cannot explain the strong decrease of the probability to observe 
a MC with moderate and large p values seen in Fig. 1. 

The axial field of the model, set to # ax i S = 1, is also well 
recovered with the fitted parameter Bq of the Lundquist field 
(Fig. 4b). Only large differences (> 20%) are obtained for glanc- 
ing encounters (p, so y/b, close to 1) or peaked magnetic field 
profiles (e.g. n = 2). 

Next, we normalize the deviation, Eq. (2), by <B>, the aver- 
age of the field strength along the simulated trajectory. Figure 4c 
shows that the deviation of the fit to synthetic data is relatively 
small unless extreme cases are considered (e.g. n = 2, see also 
Fig. 3). It is also important to notice that dev/ <B> is not a 
secure indicator of the precision of the fitted parameters: for ex- 
ample Bl fits well the synthetic observations for large p with a 
low dev/ <B> value (Figs. 3,4) while the fitted parameters are 
more biased (e.g. Bq) and/or unprecised for these glancing en- 
counters. 

Finally, Lepping & Wu (2010) found that <B>/B (p) from 
the analyzed MCs was following very well the expected relation 
from the Lundquist field except from a slight shift in ordinate 
(see their Fig. 5). This shift could be explained by n « 0.7- 
0.8 (Fig. 4d). This is an indication that the typical field in MCs 
is in between a linear force-free field and a field with constant 
axial current density. This result is in the same line of those 
obtained by Gulisano et al. (2005), where crude approximations 
were done (as assuming circular cross section, zero impact pa- 
rameter, and the orientation of the main axis given only from 
the minimum variance method). They found clues in favor of 
magnetic configurations in between linear force-free field and 
constant current models. 
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Table 1. Description of the main parameters and where they are defined. 



Synthetic models 

a bending of the boundary Sect. 4.5 

b aspect ratio of the boundary Fig. 2 

2? axis field strength on the flux rope axis Sect. 3.1 

c B z forA=l, Eq. (4) 
C = {a, b, n, r Bm in, <^ m in}> set of all characteristics 

of a model and selection parameters Eq. (8) 
n exponent defining B z and j z Eqs. (4, 5) 

x coordinate along the simulated trajectory Fig. 2 

y coordinate across the simulated trajectory Fig. 2 

z coordinate along the flux rope axis Sect. 3.1 

y/b true impact parameter Fig. 2 

Fitted Lundquist model 

a linear force-free field constant Eq. (1) 

B L Lundquist field Eq. (1) 

B estimated axial field strength Eq. (1) 

dev function of fit minimization Eq. (2) 

R flux rope radius (for B z = 0) Fig. 2 

y L estimated distance of the spacecraft trajectory Fig. 2 

to the flux rope axis 

P = yh/R, estimated impact parameter Fig. 2 

Estimated along the spacecraft trajectory 

<B> average B strength Sect. 3.2 
<B X > average B component parallel to the spacecraft Sect. 4. 1 

trajectory 

r Bx = <B X >/<B> Sect. 4.1 

co rotation angle of B across the flux rope Sect. 3.3 

Selection parameters 
^Bmin minimum average field strength to detect 

a flux rope: <B> /5 axis > r Bm i n Sect. 4.4 
co^ m minimum rotation angle of B to detect a flux rope Sect. 4.4 

Probability functions 

Pobs(p) observed probability Fig. 1 

P(p, C) theoretical probability for a model defined by C Eq. (8) 

P(p) Pip, C) integrated on a flux-rope set Eq. (9) 

P(b) probability distribution of b Eq. (9) 
^mean mean value of b Eqs. (12, 13) 
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Fig. 5. Dependence of the rotation angle, co, of the magnetic field 
component orthogonal to the axis in function of the impact parameter, 
p, found by fitting the Lundquist field to circular models (Eqs. (3,4)). 
The dots are the results obtained by Lepping & Wu (2010) for 65 MCs 
(see their Fig. 3B). 



3.3. Is there a significant selection effect with p? 

The rotation angle of the magnetic field along the simulated tra- 
jectory is weaker as the impact parameter increases (Fig. 5). 
Since an important field rotation is one key ingredient to define 
a MC, a too weak rotation angle could lead to no MC detection, 
so a bias in the probability distribution in function of p. As in 
Lepping & Wu (2010), we analyze the rotation angle of B, in 
the plane orthogonal to the flux rope axis, by taking the angle 



formed by B at each of the two boundaries (hereafter noted cS). 
As they averaged the observed B over 1 hour (their Fig. 3B), and 
to be comparable to the observations, the modeled B are aver- 
aged over 5% of the crossing length near the boundaries. For a 
typical MC duration of 20h, this implies an average over lh, so 
that our results are directly comparable to their Fig. 3B. 

For the case n = 1 the synthetic data are derived from the 
Lundquist field, so as the fitting field; it implies that the rotation 
angle has a simple expression, 2 arccos(/?), which is nearly iden- 
tical to the green curve in Fig. 5 (the differences are only due to 
the small averaging performed near the boundaries). Finally, we 
found that a broad range of n values are compatible with their 
results (Fig. 5) so that the amount of B rotation angle is not se- 
lective of different models. This result is confirmed in Sects. 4 
and 5. 

Is there a severe selection effect with the amount of rotation 
angle in observed MCs? In fact, for a rotation angle lower than 
30° no MC is observed (Fig. 5). Four MCs are observed with a 
rotation angle as low as « 50°, or lower, showing that MCs with 
low rotation angle can be detected. This value of co corresponds 
to p > 0.9 with the above models, so that a selection effect on 
MCs with low B rotation angle cannot explain the progressive 
decrease of the detection probability of MCs with p (Fig. 1). 
A way out would be to argue that significant B structures are 
frequently present within MCs, especially close to the boundary 
(i.e. for large p values), so that they can mask the lower rotation 
cases even more. However, an important rotation angle, co > 90°, 
is still present for p « 0.7 both for MCs and simulated flux ropes 
(Fig. 5). Then, it is very unlikely that the presence of B structures 
in MCs can decrease the probability of detection by a factor 3 to 
4 for p « 0.7 (Fig. 1). It implies that a selection effect on rotation 
angle cannot explain the observed probability. In the same line, 
#o remains close to # a xis (Fig- 4b) so that a selection effect on 
the field strength is expected to be weak for flux rope with nearly 
circular cross section. 

We conclude that the explored circular models cannot ex- 
plain the observations (Fig. 1). 



4. Detecting flux ropes with elongated 
cross-section 

In this section, we explore mainly the effects of the flux-rope 
boundary shape. An elliptical boundary is described by its as- 
pect ratio b (b > 1 means that the cross section is elongated or- 
thogonally to the MC trajectory, Fig. 2). We also consider bent 
cross sections in Sect. 4.5 described by an extra parameter called 
a. A rectangular cross section is also considered as an extreme 
case. These different types of cross sections allow to explore the 
space of parameters (Fig. 6) with models depending on a set of 
parameters (a, b, n). 

4.1. Expected effect of the cross-section aspect ratio 

Gulisano et al. (2007) have shown that the ratio r Bx = <B x >/<B> 
is a function of the true impact parameter y /b for a variety of cir- 
cular models (<> means averaging along the spacecraft trajec- 
tory within the MC). Demoulin & Dasso (2009) have extended 
this relationship for linear force-free models with various bound- 
ary shapes (see their Fig. 10). For an elliptical boundary, this 
relationship is summarized as r^iy/b, b,n = 1). This applies in 
particular to the Lunquist field (b = 1) and it is simply summa- 
rized as r Bx ,h(y/b) = r Bx (y/b,b= 1 , n = 1). 
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Fig. 6. Drawing defining the regions of the parameter space explored. 
The red line indicates the circular models analyzed in Sect. 3. The MC 
boundary is elliptical for the blue region and it is deformed to a bean 
shape in the green region. The two blue lines indicate the elliptical 
models analyzed in Sect. 4. Finally the purple line indicates an extreme 
case where the MC boundary is rectangular, n defines the axial electric 
current and magnetic field component (Eqs. (3,4)). A cross section 
elongated orthogonally to the spacecraft trajectory has b > 1 (Fig. 2). 
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Fig. 8. Dependence of the true impact parameter, y/b, and of the nor- 
malized deviation, dev/<B> in function of the fitted impact parameter, 
p, found by fitting the Lundquist field to elliptical boundary models. 
The parameter n describes the profile of the axial current (Eqs. (3, 4)). 



in Fig. 1). Such expectation is tested below by fitting with B L a 
variety of models with elongated cross section. 
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Fig. 7. Examples of two elliptical models (n = 0.5, 1, black dots) least 
square fitted with the Lundquist field (red curves) for a large true impact 
parameter, y/b = 0.9. B z is the axial field component, and B y is the field 
component both orthogonal to the simulated trajectory and to the flux 
rope axis. 



Next, a similar <B X >/<B> is expected when Bl is fitted to B 
(since Bl approaches the best possible B). Setting the equality 
rBx(y/b,b,n = 1) = r BXrL (p) provides a relation p(y/b,b). For a 
fixed y/b value, the derivation of this relation implies: dp/db = 
(drBx/db)/(dr B x,L/dp). Since r Bx is a decreasing function of b for 
a fixed y/b and r Bx ,L is an increasing function of p (Demoulin & 
Das so 2009), this implies that p is a decreasing function of b. 

Finally, with the magnitude of change of r Bx with b found 
in Fig. 10 of Demoulin & Dasso (2009), the value of b is ex- 
pected to strongly affect the estimated p value, so it is expected 
to strongly bias the MC probability distribution (such as shown 



4.2. Models with elongated cross-section 

We explore the space of parameters mainly with analytical mod- 
els as summarized in Fig. 6. The emphasis is set on the aspect 
ratio b since it was found to be the most important parameter 
affecting p (for a fixed y/b). We first analyze the model of Van- 
das & Romashets (2003) who derived an analytical solution of 
a linear force-free field in = 1) contained inside an elliptical 
boundary (so generalizing Bl). 

A numerical extension of the above model to cross sections 
with a bent (bean-like) shape was analyzed by Demoulin & 
Dasso (2009). They also consider the limit case of a rectangular 
cross section. It has a simple analytical expression for a linear 
force-free field (n = 1): see their Eq. (14), while their Eq. (15) 
should rather be: #r = n/2 Vl + b 2 /b. 

Finally, even if we have shown in Sect. 3 that n has a small 
effect for circular cross section, we also consider the force-free 
field with n = 0.5 and an elliptical cross section: 



(B x ,B y ,B z ) = \- 



y 



b Vl + b 2 ' VT 



4.3. Effect of the aspect ratio 



As the aspect ratio b increases, more significant deviations be- 
tween the synthetic observations of the modeled B and the fitted 
B L are present. For example, Fig. 7 shows two extreme cases 
with b = 4 and y/b = 0.9 (similar fits are obtained with lower 
y/b values). For both cases the field rotation angle, cj, is about 
120°, so large enough to be detected. However, using at typical 
value of #axis = 20 nT at 1 AU, for n = 1 such flux rope would 
not be detected, but for n = 0.5 it would since the respective 
mean magnetic field strength along the simulated trajectory is 
< 0.8 and ~ 9 nT (compare to a typical SW field « 5 nT). 

As expected in Sect. 4.1, the aspect ratio b has a strong ef- 
fect on the estimated impact parameter p (Fig. 8a,c). This ef- 
fect is much stronger than the effect of n for circular flux ropes 
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Fig. 9. Dependence of the rotation angle, oj, of the magnetic field component orthogonal to the axis in function of the impact parameter, p, found 
by fitting the Lundquist field to two models with an elliptical boundary (a,b) and one with a rectangular boundary (c). The dots are the results 
obtained by Lepping & Wu (2010) for 65 MCs (see their Fig. 3B). 
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Fig. 10. Probability distribution of the impact parameter, P(p), deduced from various models and a uniform probability distribution of y/b. The 
thin curves are without selection effect and the thick curves are for < B > /5 axis > 0.25 and a magnetic field rotation angle greater than 90°. The 
dashed black curve is a Gaussian function fitted to observations (Fig. 1). 



(Fig. 4a). The results obtained for n = 1 and a rectangular 
boundary are similar (so not shown) to the results for n = 0.5 and 
an elliptical boundary. All these results imply that p is system- 
atically biased to a lower value than the true impact parameter 
y/b, and this effect strongly increases as b is larger. 

As with circular models (Sect. 3.2), the quality of the fit of 
Bl to the synthetic data, so a low value of dev/<B>, cannot be 
used to estimate the quality of the derived fitted parameters, in 
particular of p. Indeed, even small values of dev/<B> are present 
for large p values (Fig. 8b,d) where p is the most biased com- 
pared with y/b (Fig. 8a,c). Moreover, dev/<B> has only a small 
dependence on b for n = 0.5 while p has a strong dependence on 
b. We conclude that the value of dev/<B> is not a reliable way 
to qualify the best fitting model. 

As the aspect ratio b is increased, the flux rope is stretched in 
the y direction, so one expects an increase of B y to the expense 
of B x , so an increase of the field rotation angle oj. In the models 
shown here, this enhanced rotation angle is mainly present for 
the elliptical case with n = 1 (Fig. 9b). For the two other models 
the rotation angle is almost independent of b (Fig. 9a,c). 

With a larger b value, the magnetic field can expand further 
away in the y direction, implying lower field strength (see e.g. 
Fig. 5 of Demoulin & Dasso 2009). The fit of Bl to this weaker 
field leads to a lower Bq (much lower than # a xis)- Then, we found 
that Bq is a faster decreasing function of p for larger b values. 
By contrast, <B>/Bo is found to be almost independent of b and 
more generally of the boundary shape. Then, the slightly higher 
value of <B>/Bq for MCs than for a Lundquist field, as found 
by Lepping & Wu (2010) in their Fig. 5, is mainly related to the 



B Z (A) relation, and in particular to n, Eq. (4), as found at the end 
of Sect. 3.2. 



4.4. Expected observed distribution of impact parameter 

The above bias on the estimated impact parameter, p, has impor- 
tant implications for the observed probability distribution (e.g. 
Fig. 1). More precisely, let us consider MC models with the same 
physical characteristics and observing bias (called C, which de- 
fines a set of five parameters, see Table 1). The simulated cross- 
ing is set at y/b with a distribution P(y/b). The models present 
in the interval [y/b,y/b + d(y/b)] are mapped to the interval 
[p, p + dp] with the Lundquist fit. The two probability distri- 
butions are related by: 

P(p,C)dp = P(y/b)d(y/b). (7) 

Moreover, some flux ropes could not be recognized as MCs 
because the crossing was too close to the flux-rope border. We 
include two important selection effects: a too weak field strength 
and a too low rotation angle of the magnetic field. Other selec- 
tion effects are associated to the presence of strong distortions, 
especially present when two MCs are interacting (e.g. Wang 
et al. 2003; Lugaz et al. 2005a; Dasso et al. 2009). Our models 
cannot take into account these relatively rare cases of MCs in 
interaction. For isolated MCs, the distortions close to the bound- 
ary are expected to be the strongest (weaker magnetic field and 
stronger effect of the surroundings, e.g. Lepping et al. 2007), 
therefore we select a relatively large minimum rotation angle, 
^min = 90°, while MCs are detected in observations with a min- 
imum rotation angle of « 40° (Fig. 9). The flux rope can also be 
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Fig. 11. Approximate dependence of the true impact parameter, y/b, 
in function of the estimated impact parameter, p, for bent cross-sections 
derived from Demoulin & Dasso (2009) results (derived from r Bx , see 
text in Sect. 4.5). The bending increases with the dimension-less pa- 
rameter a. Linear force-free models (n = 1, Eq. (4)) are shown for two 
aspect ratio b. The black dashed line is the relation found by fitting the 
Lundquist field to the elliptic (a = 0) model with n = 1. 



missed if its magnetic field strength is too weak. We select cases 
with <B> /#axis > ^Bmin with typically r Bm in = 0-25, since for a 
typical # ax i S value of 20 nT this implies that <B> is comparable 
to the typical field magnitude in the solar wind at 1 AU. Suppos- 
ing a uniform distribution Piy/b) (see Sect. 1), and including the 
above selection effects in Eq. (7), the probability to detect a MC 
is: 



P(P,C)-. 



d(y/b) 



dp 



(8) 



<5>/5 axis >r Bmi n & OJ>OJ 

111 



Within the studied models (Fig. 6) the weakest bias of p 
with increasing b is obtained for the linear force-free model, 
n = 1, with an elliptic cross section (Fig. 8c). It implies a 
moderate decrease of P(p,C) with p without selection effect 
(/Bmin = ^min = 0, see thin curves in Fig. 10b). The selection 
effect with oj is only present for large p values and its effect de- 
creases with b (Fig. 9b). However, this model has also the weak- 
est <B> for large p values. It implies a strong selection effect for 
TBmin = 0.25, increasing with b (thick curves in Fig. 10b). 

Increasing the axial currents (n = 0.5) or extending the cross- 
section to a rectangular shape both imply a stronger magnetic 
field for large psoa weaker selection effect. The selection effect 
with the field rotation angle remains also limited to large p val- 
ues (Fig. 9a,c). Moreover, as b is increased, the strong decrease 
of d(y/b)/dp with p further decreases the selection effects (com- 
pare thin and thick curves in Fig. 10a,c). We conclude that the 
relation piy/b) has generically a major effect on the observed 
MC distribution drawn in function of p. 

Could we interpret the observed distribution P ^(p) (Fig- 1) 
as due to oblate cross-sections? For the elliptic case with n = 1, 
none of the P(p, C) distributions, with a fixed b and selection 
criteria, is close to the observed distribution (Fig. 10b). However 
a mixture of such distributions well could be, and this will be 
analyzed in Sect. 5. For the elliptic case with n = 0.5, P(p, C) is 
very close (i. e. within the error bars) to the observed distribution 
Pobs(p) for ^ 2 (Fig. 10a), while for the rectangular boundary 
with n = 1, P(p, C) is also close to the observed distribution for 
b* 1.5 (Fig. 10c). 

4.5. Effect of a bent cross-section 

In some MHD simulations, the flux rope is strongly compressed 
in the propagation direction, such that its front region becomes 
relatively flat (e.g. Vandas et al. 2002). The cross section can 
even develop a bending of the lateral sides towards the front di- 
rection when it central (resp. lateral) parts move in a slow (resp. 



fast) solar wind (e.g. , Riley et al. 2003; Manchester et al. 2004). 
Then, Demoulin & Dasso (2009) have investigated the effects 
of bending the flux-rope boundary to a bean-like shape. This 
bending is parametrized by the dimension-less parameter called 
a. Examples of computed fields with various a values are shown 
in their Figs. 3-5. Typically \a\ needs to be larger as b increases 
to get a comparable bending. 

Does bending of the flux-rope cross section modify the prob- 
ability distribution of flux rope detection, P(p,C), versus the 
estimated parameter pi The effect of a value can be approxi- 
mately derived from the relation r# x = <B x > /<B> in function 
of y/b, a and b as summarized by the analytical expression of 
Eq. (31) of Demoulin & Dasso (2009). As in Sect. 4.1, a similar 
<B X > /<B> is expected for B and its fitted Bl field. Setting 
the equality r Bx ^(p) = r Bx (y/b,a,b) provides an estimation of 
p, named p Q (a,b,y/b), which is shown in Fig. 1 1 for a few a and 
b values. 

We also compare the estimation p Q to the result of fitting B L 
to B for a - 0, so an elliptical boundary (dashed black line). 
For b - 1, both curves are simply p - p Q = y/b (Fig. 11a), 
while for b > 1 there is a good agreement up to large y/b values 
(Fig. 1 lb). Such result could be extended to a > by applying 
the Lundquist fit to the bent models developed by Demoulin & 
Dasso (2009). Then, the analytical expression p Q (a,b,y/b) pro- 
vides an estimation of p for a broad range of {a,b,y/b} values. 
This result has a practical application: it provides a good initial 
guess of p (from observed r Bx ) for the non-linear fit of B L to B 
(so both avoiding to start in a wrong local well of dev, Eq. (2), 
and speeding up the computations). 

Figure 1 1 shows that bending of the flux-rope cross section, 
so increasing \a\, increases p for a given y/b. This is the opposite 
effect to increasing b (Fig. 8a,c). From Eq. (8), this implies a 
bias increasing the probability of flux ropes for large p, so the 
opposite to observations (Fig. 1). It is also worth noting that 
\a\ = 2 is already a very bent cross section (see Figs. 4,5 of 
Demoulin & Dasso 2009) that we expect to be rarely present in 
observed MCs. We conclude that the effect of bending the cross- 
section is expected to introduce only a weak bias to the estimated 
p value. 



5. Distribution of the cross-section aspect ratio 

In the previous section we found that the probability distribu- 
tion P(p, C), was most sensitive to the aspect ratio b. In this 
section we use this property to constrain the probability distribu- 
tion, P(b), of the aspect ratio b for the MCs observed at 1 AU. 
We end by exploring how P(b) depends on the MCs properties. 



5.1. Method 

In the following, we consider that the aspect ratio, b, is dis- 
tributed according to the probability function P(b), while the 
other parameters in C remain the same. Because of the small 
effect of a, see Sect. 4.5, we set a = 0. The expected proba- 
bility of the impact parameter P(p) is the superposition of the 
contribution of each b values according to: 



P(P) 



f 



P(b) <P(p,C)db, 



(9) 



where C = {a, b, n, rBmin, <^min) and F° P(b) db = 1, while 

-^min 

J P(p, C) dp < 1 since cases are missed with the selection on 
r Bm in and oj mm . 
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Fig. 12. Probability distributions P(p), from Eq. (9), and P(b), from Eqs. (11, 13), for the minimum of dist(rj,b mean ) as defined by Eq. (10). 
Pobs(p) of Fig- 1 is added in (a,d). Two force-free elliptical models are shown: (a-c) n = 0.5 (constant current density), (d-f) n = 1 (linear 
force-free) for r Bm i n = 0.25 and co^ m = 0. The three P(b) functions, shown with three colors, imply similar results. 



Since P(b) is contributing through an integral to the distri- 
bution P(p) in Eq. (9), and that P b S (p) (shown in Fig. 1) has 
important uncertainties due to the limited number of observed 
MCs, we can only derive a global behavior of P(b). For that, we 
limit the freedom of P(b) by selecting functions which depend on 
few parameters (b mem , b' c ), and we minimize the distance, dist, 
between P(p) and P bs(p)' 



where the coefficient in front of (Z? max - b) is computed from the 
normalisation J max P^(b, b mem ) db = 1, and 

^max 



r 



b P L (b, b mem ) db ■ 



+ 2 



(12) 



As a second possibility for P(b) we select a Gaussian distri- 
bution, limited to b > 1 : 



dist ■ 



P(b)P(p 9 c)dbj dp. (io) P Gi b,b m b> c ) = / ex p(-^^) 



(13) 



We introduce the parameter 77 in front of P(p) since P bs(p) is 
normalized with all the observed MCs (J^ P bs(p) dp = 1) while 
each P(p, C) is normalized to all the cases. Since we cannot also 
normalize P b s (p) to all cases, we let 77 as a free parameter. It is 
expected to be around 1 since the selection biases are expected 
to be small (Sect. 4). 

The generic cross-section shape of MCs is mostly unknown 
since the only shape determinations were done with a Grad- 
Shafranov reconstruction technique or by fitting the elliptical lin- 
ear force free model on a few MCs (see Sect. 1). The cross sec- 
tion has the tendency to be round (b « 1) because of the magnetic 
tension and the typically low plasma j3 found in MCs. However, 
the large pressure of the MC sheath tends to elongate the cross- 
section orthogonally to the MC mean velocity, so b > 1 . Then, 
we set a minimum value for b as b^ = 1 . Indeed, the MCs with 
b < 1 cannot be too numerous, otherwise more MCs with large 
p would be observed (because for b < 1 the bias of p(y/b) is 
reverse than for b > 1). 

We first select a simple linear function for P(b): 



Ph(b, ^mean) = 



(£max " l) 2 





(£max-£) if 1 <b<b max (11) 

otherwise, 



with 



/ 

K 



n (1 +ert{b' c )) a 
(& c -l)/(V5<r), 



yft [1 + erf (b' c )] 



+ 1 



where erf is the error function. The coefficient / was computed 
from the normalisation f Pg(^> ^mean, b' c ) db = 1. The parame- 
ter b mQm is the mean value of Pq, restricted to b > 1. The free- 
dom of Pq is expressed in function of {b mean , b' c \ rather than with 
the usual parameters of a Gaussian distribution {b c , cr} (Eq. (13)) 
in order to easily compare our results with the linear distribu- 
tion Pl (Eq. (11)). Moreover, as shown below, b mem value is 
the most stable result deduced from minimizing the function dist 
(Eq. (10)). Therefore, we set the parameter b mQm in both dis- 
tributions. For a given b mQm , the normalized parameter b' c de- 
termines the location of the maximum of Pq and the spread of 
the distribution, as follows. The probability at b = 1 divided by 
the maximal one, at b = b c , is simply exp(-Z/ 2 ) so W c describes 
how much the function Pg is peaked (b c = b' c cr V2 + 1 , then for 
b' c - 0, its maximum is at b = 1 , while it is more peaked toward 
b > 1 as b' c increases). 
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Fig. 13. Effect of r Bmin (the selection is defined by < B > /5 axis > 
'"Bmin along the simulated crossing). The coefficient 77, b mean are found 
by minimizing dist (Eq. (10)) for two force-free fields (lower curves: 
n = 0.5, upper curves: n = 1). Three probability distributions of P(b) 
are shown with different colors (for n = 0.5, the three curves are almost 
identical). 

5.2. Probability distribution of aspect ratio: P(b) 

In this section, for a given P bs(p) (Fig- 1), the function 
dist(rj 9 b mGm ), defined by Eq. (10), is minimized. We provide 
typical results for P(b). 

The function dist(rj,b m&m ) has a well defined global min- 
imum in all explored cases, see e.g. Fig. 12c,f where cuts 
through the minimum is shown in function of b mQm . For n = 0.5 
the minima are nearly at the same location (77 « 0.91 ± 0.01, 
^mean ~ 2.29 ± 0.01 for the three P(b) functions shown, while 
for n = 1 77 is larger (« 1.26 ± 0.06) and b mQm is more broadly 
distributed (from « 2.2 to 3.4). For each n value, the derived 
P(p) are all very close and fit globaly well the observations 
(Fig. 12a,d), with a comparable minimum of dist 0.036 for 
n = 0.5 and « 0.035 for n - 1). There are still some differ- 
ences: for the case n = 1, P(p) is slightly lower than P bs(p) 
for both small and large p values (p < 0.3 and p > 0.7), while 
it is the opposite for the case n = 0.5 (Fig. 12a,d). It is an in- 
dication that n is typically in between these values in MCs, in 
agreement with the result found for <B> /B Q (p) at the end of 
Sect. 3.2. 

5.3. Sensivity ofP(b) 

We compare below the results for Pl and Pq varying both the 
models (n, cross-section shape) and the selection effects (r Bm in 
and oj min ). 

The results above are derived by fitting the theoretical results 
to P bs(p) which has statistical fluctuations with the relatively 
low number (100) of MCs available. Then, we also derive the 
results from the Gaussian and linear fits (Fig. 1). The larger 
change is present for the case n = 1, and we find that b mQm 
is inside the range [2.2, 3.4] for the P(b) distributions shown in 
Fig. 12e. The range found for b mem is changed to [2.7, 3.0] when 
the Gaussian fit is used, and to [2.3, 2.5] for the linear fit. For n = 
0.5, the changes are more limited: b mean « 2.29 with P bs(p), 
changing to « 2.33 for the Gaussian fit and « 2.1 1 for the linear 
fit. We conclude that the results are weakly dependent on the 
details of the function P bs(p)- 

The selection on rotation angle, o)^, has a low effect on 
the minimum of dist(rj, b mQan ) for co m i n < 90°. The main effect 
of increasing cj m - m is to force P(p) to zero for large p values 
(Fig. 10). This effect remains in the integration on b in Eq. (9). 
For example, with oj m { n = 90°, P(p) = for p > 0.75 for both 
n = 0.5 and 1, in contradiction with Pohsip) (Fig. 12). However, 
when oj m { n is decreased to « 45°, there is only a slight decrease 
of P(p) for p > 0.9, then a> m i n around 45° is compatible with 
Pobs(p) m agreement with the minimum rotation angle detected 
in MCs (Fig. 9). 



We next explore the sensitivity of the results with rBmin se- 
lection. The elliptical linear force-free field (n = 1) is the most 
affected by changes of r Bm in threshold (Fig. 13). This is indeed 
expected from the results of Sect. 4.4, and in particular from 
what is shown in Fig. 10b. As r Bm in increases, so does the selec- 
tion effect for large p values, and lower b values are needed to 
fit the observations and a larger 77 is needed to compensate the 
selection effect (Fig. 13). At the opposite, the case n = 0.5 is 
almost independent of r Bm in since the selection affects only the 
low probability tail of P(p,C), see Fig. 10a. Similar results are 
obtained for n = 1 and a rectangular shape, with only a shift of 
^mean to « 1.56 ± 0.01 and 77 increasing a bit to 1.07, as expected 
from Fig. 10c. 

We conclude that the observed probability P bs(p) is 
mostly affected by the oblateness, b, of the flux-rope cross 
section. 

5.4. Main constrain on P(b) 

The results above are also relatively independent on the function 
P(b) selected within the explored set. In all cases close results 
are obtained from a linear and Gaussian distribution having a 
maximum at b = 1 (e.g. Fig. 12). Moreover, similar results are 
found for a Gaussian distribution more peaked around its max- 
imum, especially for the elliptical n = 0.5 and the rectangular 
n = 1 cases, i.e. changing b' c has nearly no effect on 77 and 
^mean values minimizing dist(r/, b mesm ). This is illustrated by the 
cases b' c = and 1 in Figs. 12,13. This is also true for much 
larger b' c , so more peaked Gaussian function (indeed also in the 
limit b' c — > 00, so when Pg select only b - Z? mea n)- This prop- 
erty is linked to the behavior of the functions P(p, C): the ones 
for b « Z? mean approximately fit the observations while the ones 
for larger b are too peaked to low p values and the opposite for 
lower b values (Fig. 10a,c). Then, for a distribution of b values, 
the best fit is always found around the same Z? mea n value, and the 
behavior of P(p, C) for lower b values tend to compensate those 
for higher b values. 

The above results can be modeled with the following analyt- 
ical functions: 



1 ylq o n 9 
P(p, anal.) = v ^ exp(-qp 2 b) , (14) 
erf( ^/qb) 

which approximate the behavior of P(p, C) for the n = 0.5 ellip- 
tical case with q « 1.4, and for the n = 1 rectangular case with 
q « 2. The n = 1 elliptical case has P(p, C) functions the most 
different from P(p, anal.), while still with some global similar- 
ities in their dependences with p and b. Then, and even in this 
case the results are weakly dependent on b' c (Fig. 13). We con- 
clude that the observations, summarized with Pobs(p), mainly 
determine the mean value of b, independently of the shape of 
P(b). 



6. Application to subsets of MCs 

6.1. Correlation between MC parameters 

In this section we explore the correlations between p and the 
other global parameters measured in the set of 100 MCs ob- 
served at 1 AU. In particular we find unexpected correlations. 

First, we examine the cone angle ft which is the angle be- 
tween the MC axis to the solar radial direction. The number of 
detected MCs decreases with a lower J3 angle (Fig. 14a). Still, 
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Fig. 14. Correlations of the impact parameter p with: (a) the angle between the MC axis and the radial solar direction, (b) the mean MC velocity 
(in km/s), and (c) the flux-rope radius (in AU) found with the Lundquist fit. The straight line is a linear fit to the data points (MCs). 



we find no correlation between J3 and p, showing that the cross- 
ing cases away from the MC nose (low p values) have no special 
biased impact parameter when the few cases corresponding to a 
leg crossing are filtered out (fl > 30°). This justifies the use of 
models with locally straight axis (e.g. Owens et al. 2012, and 
references therein). 

As reported by Lepping & Wu (2010), we also find no sig- 
nificant correlation between p and Bq (the deduced axial field 
strength). We agree with their interpretation that # a xis of MCs 
is expected to be spread in a large range (about a factor 10), so 
that the dispersion of # a xis is likely to mask any weak depen- 
dence Bo(p). Indeed, we find such dependence in the models. 
The dependance is weak for circular models (Fig. 4) and mod- 
erate for models with elongated cross sections. For example, 
with b = 2, #o monotonously decreases from 1 to « 0.4 for 
an elliptical linear-force free field, while this decrease is much 
weaker, only down to « 0.8, for an elliptical model with uniform 
current (not shown). Since we found two indicators in favor of 
finding typical MCs in between those models (ends of Sects. 3.3 
and 4.2), Bq(p) is expected to have a relatively weak dependence 
(from 1 to « 0.6) which can be easily masked by the large dis- 
persion of # ax i s in MCs. 

Other global parameters are not or only weakly correlated 
with p except two: V (mean velocity of the MC along the 
spacecraft trajectory) and R (flux rope radius deduced from the 
Lundquist field). The Pearson's correlation coefficient is 0.26 
and 0.35 for V and R respectively, and a linear fit also clearly 
show the trends (Fig. 15b,c). The correlation V(p) is the most 
surprising since V is measured directly from the data and it is a 
robust quantity (weakly dependent on the selected MC bound- 
aries). Such result could not be interpreted as a real velocity 
shear between the MC core and its surrounding since by its mag- 
nitude this effect would shear apart the flux rope before its arrival 
to 1 AU (and consequences of this particular behavior around 
1 AU is both not observed and not plausible). The strong cor- 
relation R(p) is also surprising. Still, we emphasis the study of 
V because R could be affected by the amount of reconnection 
achieved between the MC and the overtaken magnetic field as 
deduced by the presence of a back region in MCs (Dasso et al. 
2006, 2007; Ruffenach et al. 2012). 

6.2. Sets of MCs with different aspect ratios 

We investigate the above puzzling result by analyzing probabil- 
ity distributions, as in Fig. 1, but for MCs with a restricted inter- 
val of velocity. Due to the relatively low number of MCs, we are 
limited to a relatively coarse sampling in V. 

The probability distributions are fitted by a straight line (as 
the black line in Fig. 1) in order to decrease the statistical fluctu- 



ations inside the p bins and to summarize the distribution infor- 
mation to the slope. For an histogram of N MCs and distributed 
according to a linear function of p, the constrain that the sum of 
the probabilities is unity implies a relation between the slope of 
this distribution and the mean value of p, noted < p >, as: 

slope = 12 Ap (< p > -l/2)/(l - AT 2 ) , (15) 

with Ap being the bin size. For N slightly large (say N > 10), 
Eq. (15) shows that the slope is almost independent of N and 
simply related to < p >. It implies that the slope is a relatively 
robust quantity, even for a low number N of MCs used to build 
the distribution. The expected statistical fluctuations on < p > 
are of the order of < p > / yfN, which translate to fluctuations of 
the slope « (slope + 6Ap)/ V^V « 0.1 for a slope « -0.2 (Fig. 1), 
Ap = 0.1 and N = 16. 

Next, we ordered the MC data according to growing values 
of V, and computed the evolution of the slope for N MCs pro- 
gressively shifting to higher V values. With N = 16, fluctu- 
ations of the slope are < 0.1, as expected. There is a sudden 
change for V above ^ 550 km/s (Fig. 15a). A similar result is 
obtained for larger N values, with less fluctuations, but with a re- 
duced dynamic (in both axis directions). Indeed, separating the 
MCs in two groups shows two different distribution functions 
(Fig. 15b,c). Similar results are found when the above ordering 
with V is replaced by one with R. 

With the results of Sects. 5.2 and 5.4, we interpret this result 
as the presence of two main groups of MCs. The slower ones, 

V < 550 km/s, which are also the most numerous (84 MCs), 
have an oblate cross section with a mean aspect ratio between 2 
and 3 depending the model used, similar to the full set of MCs. 
However, the faster MCs at 1 AU have a nearly flat distribution, 
so they are mostly round whatever model is selected (within the 
explored ones). It would be worth checking this conclusion with 
more MCs since this group is limited to 16 MCs. These MCs 
are also typically larger and with stronger magnetic field since 

V has a correlation coefficient of 0.32 with R and 0.68 with Bq 
(for the full set of 100 MCs). Indeed, a variation of the slope of 
*Pobs(p) with MCs ordered with R was found to be similar than 
with V (Fig. 15a). This is not the case with Bq since there is no 
significant correlation between Bq and p (Sect. 6.1). 

Why faster and larger MCs would have typically nearly 
round cross section? At first thought, a faster MC would imply 
a larger snowplow effect, plausibly generating a larger sheath 
which can compress more the flux rope, inducing a flatter cross 
section. However, the velocity is measured at 1 AU and the 
above result could rather mean that those faster MCs were in 
average less decelerated than others, so that the distortion from 
the surrounding solar wind was less important than for other 
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Fig. 15. Properties of impact parameter distributions P bs(p) for dif- 
ferent MC groups, (a) Slope of the linear fit of P b S (p) when MCs are 
first ordered with a growing mean velocity V, then binned in groups of 
16 cases (running grouping with increasing V). The three curves rep- 
resent this slope, with the black line corresponding to the mean V of 
each group, and the blue (resp. red) line corresponding to the minimum 
(resp. maximum) value of each group. The horizontal dashed line is the 
slope for all MCs (slope of the black line in Fig. 1). The horizontal axis 
is scaled with the logarithm of V. (b,c) Probability distribution, Pob S (p), 
as in Fig. 1, with MCs separated in two groups according to their mean 
velocity V. 

MCs. Moreover, faster MCs spend less time from solar erup- 
tion to their arrival to the point where they are in situ observed, 
and then distortion mechanisms are expected to be less effective 
to operate. Another, and plausibly complementary answer is that 
the faster MCs have typically a stronger magnetic field, so that 
the magnetic tension is stronger and keep the cross section more 
round. 



7. Conclusions 

The MCs observed at 1 AU are classically fitted with a Lundquist 
model (Lepping et al. 1990). In the set of 120 MCs analyzed, 
only 11% (13/120) of the MCs could not be satisfactory fitted 



(either the flux-rope handedness could not be determined, either 
the fit was not converging), while 6% (7/120) of the MCs are 
crossed to far away from the nose to provide reliable fit results. 
For the remaining 100 MCs the fit provides an estimation of 
the impact parameter (p). The observed probability distribution, 
Pobs(p), of these MCs is found to decrease strongly with p (Lep- 
ping & Wu 2010, and Fig. 1). Compared to an expected almost 
uniform distribution, this could imply that about half of the MCs 
are not detected by in situ observations. Is this decrease due to a 
strong selection effect, like on the magnetic field strength and/or 
the amount of field rotation angle? Or perhaps are the MCs ob- 
served only in about one third of ICMEs because more criteria 
are used to define ICMEs than MCs? Moreover, several of the 
less restrictive criteria used to identify ICMEs are expected to be 
independent of the impact parameter (such as temperature, com- 
position and ionization level). In order to answer these ques- 
tions we explore the parameter space of flux rope models with 
force-free fields. We simulate spacecraft crossings and perform 
a least- square fit of the synthetic data with a Lundquist field, us- 
ing the same procedure as for observations of real MCs. The fit 
provides an estimated impact parameter p that we compare to 
the true one known from the synthetic model. 

For models with circular cross sections, we found that selec- 
tion effects with magnetic field strength and field rotation angle 
are present only for large p values, so they cannot explain the 
gradual decrease of P bs(p)- This result is found for a broad va- 
riety of magnetic field profiles ranging from nearly uniform to 
peaked field strength across the flux rope. 

Next, exploring non-circular cross sections, we found that 
the aspect ratio, b, of the cross section is the main parameter af- 
fecting the estimated impact parameter p. For flux ropes flatter in 
the propagating direction (corresponding to b > 1), p is more bi- 
ased to lower values, compared to the true one, as b is increased. 
This effect implies simulated distributions P(p) which are close 
to observed ones with b « 2 for an elliptical model with uniform 
axial current density. For linear force-free fields with elliptical 
cross-section, p is less affected by b, but the field strength de- 
creases more rapidly away from the flux-rope axis, so that the 
selection effect on the field strength enhances the dependence of 
P(p) on b. 

We also explore other effects which can bias the probability 
distribution of p. We found that bending the cross section in a 
bean-like shape has a small effect on the estimated p. A much 
larger effect is present if the cross section is set broader than an 
ellipse at large distance from the axis. An extreme case is a rect- 
angular cross section. In that case, the linear force free model 
corresponds to an even more biased p than the above elliptical 
model with uniform axial current density, and b « 1.5 is suffi- 
cient to reproduce the observed distribution P bs(p)- Finally, we 
found that for all the models explored, the rotation angle along 
the spacecraft trajectory is above 90° except for large p values (at 
least p > 0.7). Then, a selection effect on this parameter cannot 
explain P \>s(p). Furthermore, only a selection criterium around 
40° can lead to a computed P(p) in agreement with P h S (p) for 
large p values. This is in agreement with the minimum rota- 
tion angle found in the set of MCs analyzed by Lepping & Wu 
(2010). 

We conclude that the observed distribution P bs(p) is mainly 
shaped by the oblatness of the MC cross section, with some con- 
tribution of a field strength selection when the flux rope is close 
to a linear force-free field. Still, even in this last case, typically 
more than 70% of the flux ropes are expected to be detected. 
Even adding unfrequent cases not detected because of very large 
perturbations (so that the field rotation is not detected), or with 



Article number, page 13 of 15 



a crossing within a leg, or MCs strongly in interactions, this im- 
plies a low amount of undetected flux rope, well below 2/3. So 
we conclude that a large majority of flux ropes are expected to 
be detected. The non-MC ICMEs are either encountered outside 
the flux rope limits or they would contain none. 

We also get results beyond the initial questions. The main 
dependence of P(p) on the aspect ratio b allows to constrain 
a key property of the distribution P(b) for MCs: the mean of 
the aspect ratio. With an elliptical model with uniform current 
density, P bs(p) sets the mean of b nearby 2.3 independently of 
the broadness of the distribution. This last property is approxi- 
mately kept for a linear-force field, but the mean of b is shifted to 
around 3 with a slight dependence on the amount of the selection 
effect of the field strength. Then, we conclude that the observed 
'Pobs(p) implies that MCs are moderately oblate at 1 AU, at least 
in average. 

We further analyze the observed MCs by separating them in 
groups with different physical parameters. In contrast to most 
MCs, the faster MCs (above « 550 km/s) have a flat P bs(p) dis- 
tribution. It implies that the faster MCs, which typically have 
also both larger radius and field strength, are nearly round, while 
the slower ones have typically the above mean oblateness. Fi- 
nally, we find two results indicating that the typical magnetic 
field profile in MCs is in between a linear-force-free field and 
one with a constant axial current density: 

- First, the mean field strength observed along the spacecraft 
trajectory is systematically above what is predicted by a linear 
force-free field, but below the prediction given by a constant 
current model, and this is independent of the aspect ratio of the 
cross section, in agreement with a previous study (Gulisano et al. 
2005). 

- Second, the distribution P(p) computed with a distribution 
of b, derived to fit *P bs(p), shows systematic biases, both at low 
and large p values, with opposite tendency for both type of mag- 
netic fields. This conclusion is also coherent with the flatter field 
strength profile found in MCs compared to a Lundquist field. 
Then, both the current distribution and the oblateness of the flux 
ropes contribute to a relatively flat profile of the field strength. 
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